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ABSTRACT 

We propose that high frequency quasi-periodic oscillations (HFQPOs) can be 
produced from randomly-formed X-ray bursts (flashes) by plasma interior to the 
ergosphere of a rapidly-rotating black hole. We show by direct computation of 
their orbits that the photons comprising the observed X-ray light curves, if due to 
a multitude of such flashes, are affected significantly by the black hole's dragging 
of inertial frames; the photons of each such burst arrive to an observer at infinity 
in multiple (double or triple), distinct "bunches" separated by a roughly constant 
time lag of Atiag/M ~ 14, regardless of the bursts' azimuthal position. We argue 
that every other such "bunch" represents photons that follow trajectories with 
an additional orbit around the black hole at the photon circular orbit radius (a 
photon "echo"). The presence of this constant lag in the response function of 
the system leads to a QPO feature in its power density spectra, even though the 
corresponding light curve consists of a totally stochastic signal. This effect is 
by and large due to the black hole spin and is shown to gradually diminish as 
the spin parameter a decreases or the radial position of the burst moves outside 
the static limit surface (ergosphere). Our calculations indicate that for a black 
hole with Kerr parameter of a/M = 0.99 and mass of M = IOMq the QPO is 
expected at a frequency of z/qpq ~ 1.3 — 1.4 kHz. We discuss the plausibility and 
observational implications of our model/results as well as its limitations. 



Subject headings: accretion, accretion disks — black hole physics — X-rays: 
galaxies — stars: oscillations 
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Introduction 



Following the inspiring observational discoveries of quasi-periodic oscillations (QPOs) 
from compact objects [see, e.g., van der Klis 2000 for neutron star low-mass X-ray binaries; 
Strohmayer 2001a,b and Cui et al. 1999 for stellar-mass black hole systems; Strohmayer 
et al. 2007 for ultraluminous X-ray sources (ULXs)], a number of theoretical scenarios have 
been proposed to explain the physics behind the QPO observations. These models are based, 
among others, also on the dynamics of relativistic accretion disks, especially for the QPOs 
of systems thought to harbor black holes. For example, it has been proposed that the X-ray 
modulation giving rise to the observed QPOs can be produced at the precession frequency of 
accret ion disks due to relativi stic dragging of inertial frames around rapidly-rotating black 
holes (ILense fc Thirring Il918f). with the QPO frequency used t o estimate the range of the 



hole' s spin parameter (e.g. ICui et al.l 1 19981 ; lAschenbachl I2004J ; ISchnittman fc Bertschinger 
20041 : ISchnittman et al.l [20061, for XTE J1550-564, GRO J1655-40, GRS 1915+105, Cyg 
X-1, and GS 1124-68). On t he other hand, resonance freq uencies among various disko- 
seismic oscillation modes (e.g. iNowak et al.l 119971 ; lKatdl200ll ) have been invoked to explain 
the observed 2:3 frequency co r nmensurability in XTE J1550-564 and GRO J1655-40 (e.g 



Abramowicz &: KluzniakI 1200 ll : iKluzniak fc Abramowica 1200 ll : ISchnittman fc Bertschinger 



20041 : bonmeJ I2OO7I ) More recently, lAoki et al.l (|2004l ) argued that radially propagating 



shocks in relativistic accretion disks could produce the observed QPOs from black hole sys- 
tems (GRS 1915+105 and GRO J1655-40). Finally, considerations of inhomogeneities of 
accretion disks (e.g., due to local magnetic flares a nd /or orbiti ng clumps) have led some 
authors to investigate orbiting hot spot models (e.g. lKaraslll999l ). The models in the above 
(necessarily incomplete) list can provide frequencies in agreement with those of the observed 
QPOs by the judicious choice of some of the systems' dynamical parameters. The QPO 
features in these models rely on an underlying oscillatory behavior of the light curve which, 
however, can be easily obscured in the presence of noise, even though the QPO is clearly 
present in the power spectra. 

In this paper, we discuss a process for QPO formation that relies not on an underlying 
(but noise-obscured) oscillation, but on an "echo" of the input signal (the QPO behavior 
then follows from a well known theorem of Fourier analysis). As such, we show that QPO 
features are possible even for a totally random signal, which in itself differs little from white 
noise. We show that such an "echo" (and the concomitant QPO) is possible for photons 
emitted by an accretion disk whose Innermost Stable Circular Orbit (ISCO) reaches within 
the ergosphere of a rotating black hole, even if the photon emission is random in time, orbital 
phase and isotropic in the plasma frame. Our analysis of this process indicates that these 
features imply the presence of, unobserved as yet, high frequency QPOs (HFQPOs), which 
rely crucially on the dragging of inertial frames and are absent for slowly rotating black 
holes. 
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Our paper is structured as follows: In §2 we provide a general description of the model, 
details of the photon kinematics and the response function of the system. In §3 we give a 
prescription for constructing stochastic model light curves and show that their power density 
spectra (PDS) exhibit the QPO features as anticipated. Finally, in §4 we review our results, 
make contact with observations and discuss prospects of future work. We show the derivation 
of the radial null geodesic equation in the Appendix. 



2. Description of the Model 



The model we consider consists of the standard geometrica lly thin, optically thick ac 



cretion disk (e.g. iNovikov Sz Thornelll973l : iPage fc Thornelll974l ) that extends to the ISCO. 



We assume the production of instantaneous X-rays in short bursts (short compared to the 



connection of magnetic field loops anchored in the disk (e.g. 


Galeev et al. 


1979; 


Haardt et al. 


1994; 


Poutanen & Fabian 


1999: 


Navakshin & Kazanas 


2001; 


Czernv & Goosmann 


2004 


), or 
sider 


of standing shocks (e.g. 


i^ukumura et al. 


2007 


). To simplify our calculations we com 



observers at small latitudes and we can thus restrict with good accuracy the computation 
of photon orbits on the equatorial plane {6 = 7r/2). We assume that the X-ray flares occur 
along the circumference of a radially thin annulus randomly both in position and in time. 
The small height of the disk and the source of the X-ray flares guarantee that except for the 
photons intercepted by the disk, the remainder can reach unimpeded the observer at infinity 
at nearly equatorial orbits. We assume that the photons are emitted isotropically in the ro - 



tating plasma rest frame [this entails some subtleties (see e.g. iFukumura fc Kazanas 1120071 )] 
and their trajectories are influenced by both the motion of the emitter and the dragging of 
inertial frames (the marginal orbit is inside the ergosphere for sufficiently large values of the 
black hole spin). We collect the photons at a large radial distance {voo/M = 600 where M 
is black hole mass) as a function of time for different relative positions between the X-ray 
flare source and the observer to compile the response function of the system for observers at 
small disk latitudes. 

To compute the photon orbits we adopt the Kerr metric in Boyer-Lindquist (BL) co- 
ordinates (r, 6, 0) and geometrized units {G = c = 1 where G is the gravitational constant 
and c is the speed of light). Distance and time are normalized by black hole mass M [i.e. 
r = 1.5 X 10^ {M/Mq) cm and t = 5.0 x 10"^ {M/Mq) s]. The accretion disk extends from 



an unspecifled outer radius to an inner radius (of ISCO) at r 



within which the matter 



freely spirals in toward the event horizon. Hence, we assume that the X-ray source at (r^, 0^) 
in the disk region follows the Keplerian motion, while the source inside the ISCO must have 
radial and azimuthal motion with its energy and angular momentum at the ISCO being 
conserved subsequently. Plunging motion measured in a locally non-rotating reference frame 
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(LNRF) is described by 



A. 



A. 



(2) 



where = — 2Mrs + a^, = (r^ + a^)^ — a^A^ and a denotes a dimensionless Kerr 
parameter. The angular velocity of the source is given by VLg = u'^/u^ while u describes the 
angular velocity due to the dragging of inertial frames, and the four-velocity components, 
M*, and u"^, satisfy the normalization condition {u^u^ = —1). For Keplerian motion we 
have ^Ig = = M^/^/(aM^/^ + Vg^'^) and u"^ = 0. We numerically solve the exact null 
trajectories by specifying the position of a burst (r^, 0s) and the photon's impact parameter 
b (or specific axi al angular momenturn ), employing the following equations of motion (see, 
e.g.. Appendix in IChandrasekharl Il983l ) 



dt 
dX 



1 / . , 2a^M 

— + + 

A V r 



2aMb\ 



d\ 
dr 
dX 



1 

A 



2aM 



1 - 



2M 



^1/2 cos^ 6 (±rA| sec 6\ - 2aMA^/^ tan 5) 
^1/2 [^3(^ _ 2M) + a2(r2 - 2M^) + 2a'^M^ cos(25)] 



(3) 
(4) 
(5) 



2Mr + a2, ^ 



(r + 2M), A is the null affine parameter, and S 



where A = r 

measures photon's emission angle between propagation direction and radial direction (e.g. 
Misner et al.lll973l . p. 675) in the LNRF (as opposed to the rotatin g fluid frame). See Ap- 
pend ix for the derivation of radial component of null geodesies (also iFukumura &: Kazanas 
20071 . for a similar discussion of poloidal geodesies). We must take into account the azimuthal 
beaming effects (photon focusing) on the emitting angle 6 at the source position (r = r^), 
particularly in the innermost regions (r > r^) where the Doppler effects become significant. 
This is defined as 



6 = arctan 



sm a 



7(f + cos a) 



(6) 



where a measures photon's emission angle between propagation direction and radial direction 
in the rest-frame of the source (as opposed to LNRF), v denotes the total three- velocity of 
the source measured in LNRF, and 7 = 1/a/1 — v'^. We stress that both quantities, a 
and 5, are defined as photon's angles between emission vector and radial component vector 
in two different frames, which are related through the Doppler boost by equation (jS]). In 
our computations a (rather than 6) is chosen uniformly in fluid frame for isotropic photon 
emission. In the Appendix we explicitly show the definition of 6 and its relation to the 
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Fig. 1. — Sample signals from two random X-ray bursts (out of a total of A*";, = 6000) at (a) 
(j)s,i = 0° and (b) = 180° for a/M = 0.99 and r, = r^s. In each phase A^ph = 200,000 
photons are isotropically emitted. Arrows indicate a roughly constant time lag Atiag/M ~ 14. 
Relative positions of the observer, source and the black hole are illustrated. 



photon's impact parameter b. Note that 5 a as v ^ (i.e. no beaming case) as expected. 
We then consider locally isotropic emission from each burst in the fluid frame and trace the 
individual photon ray by solving the geodesic equations in Kerr geometry. 

For a given source position (r^, (ps,i) (where = for the source located between the 
observer and the black hole) we collect the times tj of the photons arriving at r = r^o = 
600M within an angular bin of A0 = 1° centered at the observer, to produce the response 
function of the system for the specific source location Ii{t]rs,(ps,i) [photon/time]. A total 
of A^ph = 200, 000 photon rays were computed from a single source. By varying the phase 
(f)s,i of the source relative to the observer between and 27r we then compiled the response 
function of the system for all 0^ j for a given r^, and the obtained photons were binned in 
time bins of size M. Samples of this response for (a) (j)s,i = 0° and (b) 180° are given in 
Figured] for a black hole with a/M = 0.99 and emission from the marginally stable orbit, 
i.e. Ts = ~ 1.45 M. One can plainly see that the response function consists of a series of 
pulses separated by approximately Aiiag/M ~ 14, regardless of the azimuthal position of the 
source. These represent photons arriving at the observer either directly or correspondingly 
after one or more orbits around the black hole. Also, although the observer is physically 
closer to the source in (a) than in (b), he/she detects the signals emitted in (b) earlier. This 
is due to the combination of the strong Doppler beaming due to the rotation of the source 
and the dragging of inertial frames that send most photons in case (a) around the black hole 
rather than directly into the observer's line of sight. Sample of photon trajectories emitted 
from a source with (j)s,i = 0° (corresponding to Fig. [TK) reaching the observer are shown in 
Figure [2l 

In fact. Figure [2] encapsulates the essence of the effect discussed herein: Because of the 
frame-dragging effects for a source within the ergosphere, even those photons with negative 
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angular momenta (going backward in the local frame with respect to the hole's rotation) are 
forced to propagate in the same sense as the rotating hole (the photon escaping directly to 
the observer in the upper right quadrant of Fig. [2] is one of them). All other photons can 
reach the observer (at phase (ps^i = 0° in Fig. [2]) only by moving in the direction of the hole's 
rotation; as such they can do so by moving by an angle that is a fraction g of 27r (g = 
for dashed ray; q = 1 for gray ray, and g = 2 for solid dark ray) to produce the response 
of Figure [1^. As we show explicitly in the next section, it is the constancy of this time-lag 
that is responsible for the presence of QPOs in the system. Such a constant lag is almost 
absent, however, for photon sources outside the ergosphere or in the Schwarzschild geometry, 
leading to qualitatively different results in these cases. 

3 
2 

1 


A 

-2. 
-3 

-3-2-10123 

Fig. 2. — Representative photon rays corresponding to the response function in Figure [1^. 
The trajectory in the upper right quadrant (dashed curve) is a backward emitted photon. 
The gray and dark curves (solid curves) represent photons that reach the observer after 27r 
and An radians, respectively. 

To be sure, even in the Schwarzschild geometry one expects that some photons can reach 
the observer after going around the hole by an angle that is an integer-multiple of 27r, but their 
number is much too small to make an observable contribution to the response function. In 
support of this claim we present similar computations in a Schwarzschild geometry in which 
the frame-dragging effect is absent and the marginally stable orbit is larger {r^s/M = 6). The 
response functions for (a) (j)s,i = 90° and (b) 165° are depicted in Figure [3]for A^ph = 200, 000 
photons. It is apparent that the response function in this case comprises a very large peak 
at time t/M ~ 610, a much smaller one at t/M ~ 632 and two additional even smaller 
ones separated from this last one by Atiag/M ~ 16 and Atiag/M ~ 32 respectively, i.e 
indicating the presence of a constant lag in this particular phase too. In order to get a better 
understanding of the orbits responsible for these peaks and their relation to the corresponding 
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Fig. 3. — Same as Figure [T] but for a/M = and = = 6M. No constant time-lag is 
seen. 



orbits of the Kerr geometry we did obtain and plot in Figure H] the orbits of the photons 
corresponding to these peaks. Thus the largest peak corresponds to photons that arrive 
at the observer directly from the source (in the retrograde direction with respect to the 
source direction), while the smaller one from photons that go around the hole by an angle 
A(f) ~ 37r/2 in the prograde direction. The next much smaller peak at t/M = 632 + 16 is due 
to photons that arrive at the observer after going around the hole by an angle A0 ~ 57r/2 
in the retrograde direction and the one at t/M = 632 + 32 after an angle A0 ~ 77r/2 in 
the prograde direction. As the azimuthal angle (f)s,i approaches 180°, both the two largest 
and two smallest peaks move closer (as expected by symmetry) to obtain approximately the 
same normalization and merge for 0^ ^ ~ 180°, indicating that their time-lag is in fact not 
constant but depends strongly on the value of the phase angle (l)s,i, unlike the situation for 
the Kerr geometry of Figure [H where the lags are independent of the source phase. This is 
a crucial point in our model. 



3. Model X-ray Light Curves and Timing Analysis 



We have used the response function constructed as above to produce synthetic light 
curves. In order to avoid the introduction of QPO features by the rotation of the emitting 
sources, we have made the flare emission instantaneous and positioned the sources at random 
values of (f)s,i, uniformly between and 27r; we also introduced the flares randomly in time 
using the following prescription of the time interval ATj between the i-th and {i + l)-th 
bursts 

ATi = fx |ln{rnd(0, 1)}| , (7) 

where T is a mean timescale between bursts and rnd(0, 1) denotes a random number between 
and 1. We have chosen T = fT^rh where Torb('"s) = 27r(rs^^ + aM^/^)/M-^/^ is the Keplerian 
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Fig. 4. — Representative photon rays corresponding to the response function in Figure (Sfe- 
The trajectories in the upper right quadrant (dashed curves) are backward emitted photon. 
The gray and dark curves (sohd curves) represent photons that reach the observer after 
77 + 7r/2 and Svr + 7r/2 radians, respectively. 

orbital period at r = and / (/ ~ 1) characterizes the frequency of burst occurences relative 
to the orbital period Torb (the precise value of / is in fact not relevant for the shape of the 
power spectrum we are interested in). 

Each burst produces a characteristic signal Ii{t] r^, (pg^i) unique to that specific azimuthal 
position (j)s,i with /j(t; r^, (j)s,i) denoting the intensity of the signal of the i-th burst as shown 
in Figures [T] and [31 The entire bolometric light curve then is the incoherent superposition of 
similar signals corresponding to the randomly produced flares, i.e. 



where Nf, (set to be 6,000) is the total number of X-ray bursts used in a given model light 
curve. Note that each flash contains A^ph = 200, 000 photons. 

Using the above prescription we have produced sets of model light curves varying the 
parameters Nh and / for the same values of the black hole parameters used to produce 
Figures [1] and [31 The simulated light curve seen by the observer is shown in Figure [51 where 
Nf, = 1000 bursts are considered in (a) and the first 100 bursts are extracted from (a) in 
(b). They exhibit no apparent periodicity, given the random positions (j)s,i and times of the 
induced flares ATj. We find that the characteristic appearance of the light curve essentially 
remains the same for the Schwarzschild cases too. 




8 
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We have then used standard analysis tools for the Fourier transform F{uj) to produce 
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(a) 1000 bursts (b) 100 bursts 




t/M(x100) t/M 

Fig. 5. — Two sample light curves generated as discussed in the text, (a) The light curve of 
a series of Nb = 1000 bursts from sources at radius r = = 1.45M for a black hole with 
a/M = 0.99 and / = 1. (b) The light curve of only the first 100 bursts of the same series. 
It is apparent that these are totally random. 

the PDS 1^(07)1^ of these light curves along with the Autocorrelation Function (ACF) R{t). 
These are shown in Figure [6] for various source radii, i.e. > rms, = and a plunging 
orbit with < r^^g, for a black hole with a/M = 0.99. A characteristic peak is clearly seen 
in the PDS at angular frequency of cuqpoM ~ 0.4 and its higher multiple modes, which for 
M = IOM0 corresponds to a peak frequency z/qpo ~ 1.3 — 1.4 kHz. The ACF exhibits also 
a prominent peak at the equivalent time lag of r/M ~ 14 — 15 and perhaps encapsulates 
better the underlying physics behind the QPO, i.e. the presence of a well-defined time lag 
(i.e. an "echo") in the response of this system. This is best seen, e.g., in Figures EH-d. 
The QPO pattern (peak frequency and width) remains essentially the same regardless of the 
source position, as long as it is within the ergosphere. On the other hand, as approaches 
the static limit (r = rgtatic = 2M near the equator), the QPO signal begins to disappear (see 
Figs. [6^ and f), since the time-lags between the peaks of the response function is no longer 
constant, as discussed earlier. 

It should be noted that the choice of the times of flare injection [equation ([7])] guarantees 
that in the absence of the lags discussed above the corresponding PDS would be just white 
noise. However, because our timing resolution in this numerical experiment is M, we have 
approximated each "bunch" of photons received at a given angle and a given resolution 
interval as a Gaussian of area normalized to the number of photons received in each interval 
and FWHM equal to M. It is precisely this finite width of our "shots" that leads to the 
high frequency cut-off of the PDS given in Figures [6|i-d. Of course, at the low frequencies 
the PDS is that of white noise. 

In order to study quantitatively the QPO feature in this model we have also examined 
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its dependence on the black hole rotation a for = r-^is in Figure [3 To do so we select two 
values for the black hole spin, namely a/M = 0.995 and a/M = 0.9428. The first was chosen 
to verify that the effect becomes more prominent with the increase of the black hole spin; the 
second value was chosen by the requirement that the equatorial boundary of the ergosphere 
be equal to that of the ISCO, i.e. that = Tstatic = 2M. The first case yields clearly 
prominent QPO and peaks at the appropriate lags in the ACF. However, these features 
become almost statistically insignificant as the spin parameter approaches the critical value 
of a/M = 0.9428 (for which = r^s = '"static) for the reasons discussed in the previous 
paragraphs of this section. 

We have also performed a similar calculation for a Schwarzschild black hole case for 
fs = fms = 6M to examine whether the QPO features we have found in Kerr geometry may 
also be produced in the absence of black hole rotation. In this case the PDS (not shown 
here) exhibits no QPO features similar to those of the Kerr case (as expected). The reason 
is the absence of a constant time-lag between the two major peaks of the response function 
(see Fig. [3]); the normalization of the two minor peaks is just too small to make a difference 
in the PDS, despite the presence of roughly constant lags between them at some phases. The 
corresponding ACF also showed no particular coherent timescale. We have repeated similar 
computations for various source radii (r^ > and < rms) for Schwarzschild case and 
found that in no case the QPO features seen above were produced. As discussed above this 
is due to the absence of a constant time lag in the system response. We will come back to 
discuss this issue in §4. Of course, should there be a reason that a particular phase is favored 
as the site of the X-ray flares considered by our model, then one would expect QPOs to be 
present for sources in a Schwarzschild geometry too. 

Finally, we have checked that the feature obtained is not simply an artifact due to 
photon statistics by changing the number of bursts and the random number sequence used 
and indeed nearly identical results were consistently found. Note that for a given black hole 
rotation a, there are no free parameters other than the radial position of the flares r^, and 
the mean timescale of each burst T, neither of which changes the value of the lag in the 
response function. 



4. Discussion 

We have presented above a process that can lead to QPO features in the PDS of accreting 
black hole systems. This process is to a certain extent different from the typical QPO models 
in the literature in that it relies on the "echo" of a (however random) signal for the emergence 
of the QPO in the observed PDS. By contrast, most QPO models invoke the presence of 
an underlying (but perhaps obscured due to noise) oscillation in the emitted photon flux. 
Whether this qualifies it as an alternative QPO model or not is, in our view, a matter of 
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convention. However, it provides a mind-set concerning the QPO phenomenon that deviates 
from that of the more conventional considerations which may lead to novel future insights 
on the subject. The difference of these two distinct QPO notions is perhaps exemplified best 
mathematically by the corresponding ACF, which in one case has an oscillatory behavior 
while in the other the double peak form presented in Figures [6] and [71 

As discussed in the previous section, the "echo" of the proposed model is the result 
of the black hole angular momentum a and the ensuing dragging of inertial frames on the 
trajectories of the photons emitted within its ergosphere; as such, this "echo" is absent in 
black holes of sufficiently small values of a. This is not the first time that frame- dragging 
has been invoked to account for a certain aspect of the QPO phenomenon. For example, 
the same effect and the accompanying Lense-Thirring disk precession has been invoked to 
account for the dependence of the intrafrequency correlations of sev eral low frequency QPOs 
( IStella fc Vietri I Il998l : Ivan der Klid l2000l : ISchnittman et al.l l2006l ): the effect we describe 
herein is different in that it affects the orbits of individual photons, rather the orientation of 
the entire disk. 

One should further note that, because the "echo" lags that produde these QPOs depend 
only on the background geometry of the accreting black hole (they are roughly equal to 
the length of the circular photon orbit for this geometry), they (and also the resulting 
QPO frequencies) are independent of the source flux (the accretion rate). On the other 
hand, given that the QPOs of accreting neutron stars and most of the O. POs of accreting 
black holes do depend, in general, on the source flux ( van der KlisI 2000l ). they are likely 
due processes different from that described above. Nonetheless, as discussed in §1, QPOs 
that do not dep end on the source f lux have been discovered in a number of galactic black 
hole candidates f Strohmayer 2001al ,bl). However, the frequencies of these latter QPOs are 
much too small (fqpo ~ 100 — 400 Hz) to be attributed to the process described herein 
[z/QPO — 1400 (IOMq/M)], unless our understanding of strong gravity physics is in signiflcant 
error. 

Up to this point, the discussion of our model considered flares whose duration is much 
shorter than the orbital period of the accretion disk near its ISCO. Clearly, the "echo" 
process discussed herein is present for longer flare durations, even if the latter exceeds the 
local orbital period. Howev er, in this case the model becorn es very similar to those of orbiting 
hot spots discussed earlier (lKaraslll999l : ISchnittman II2005I ). We have produced a number of 
model light curves assuming the flare duration to be several times the ISCO orbital period 
(i.e. just like in the orbiting hot spot model) and computed the corresponding PDS and 
ACF. In this case, the geometrical echo, while still present, makes only a small contribution 
to the PDS, whose QPOs are now dominated by the orbiting hot spot periodic motion. 

While the process we discussed above is quite robust (it depends only on the geometry 
of the accreting object), its potential observability depends on a number of factors. To begin 
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with, our calculations were performed exclusively on the equatorial plane of a Kerr black 
hole; therefore, our present results are valid only for disks that are geometrically thin and for 
observers at relatively low latitudes. For observers at moderately higher latitudes (or thick 
disks), our orbit calculations have to be supplemented with those for the ^—coordinate (to 
follow the poloidal motion) that have been omitted so far in this work. However, the essence 
of the effect we discussed above, i.e. forward thrust of all photon trajectories produced 
within the ergosphere, whether equatorial or not, is expected to be present in that case 
too and thus qualitatively we expect the same phenomenon to be conditionally observable 
for non-equatorial source/observer configurations. Our preliminary calculations indicate the 
presence of such QPO for observers at latitudes at least as high as ~ 30° from the disk 
mid-plane. Quantitative analysis of this aspect of the problem, which bears the application 
of these ideas to thick disks, will be presented in a future publication. 

Another limitation of the model discussed in the earlier sections is that of the source 
radius r^. The response frames of Figure [1] were computed assuming the source to be on the 
ISCO at Tg/M ~ 1.45. Changing the source radius leads to a quantitatively similar behavior, 
i.e. the ACF still exhibits a peak at the same value of the lag r/M ~ 14, as long as the 
source of photons is located within the static limit (ergosphere), i.e. for r^/M < 2 near the 
equatorial plane. Interestingly, for the value a/M = 0.99 used in these calculations the ACF 
peak at r/M ~ 14 disappears gradually as the source radius increases past r^/M > 2, i.e. 
as the source moves outside the ergosphere, a fact which suggests that the effect we have 
presented is in fact due to the strong dragging of inertial frames. Equivalently, we anticipate 
the gradual disappearance of this peak in the ACF as the value of the black hole spin a 
decreases down to a critical value of a/M ~ 0.9428 (where the radius of ISCO becomes 
equal to that of static limit) and the ergosphere shrinks to leave much of the disk outside. 
While the frame-dragging is still operative even at radii outside the static limit, it is simply 
not effective enough to cause significant azimuthal beaming. 

We have in fact checked that the effect we consider is not due to the source proximity 
to the horizon (a situation possible without free fall onto the hole only if the latter is rapidly 
rotating) by computing the orbits of photons from radii as small as r^/M ^ 2.1 by matter 
in-falling (in the plunging region of < r.^^) onto a Schwarzschild black hole. In computing 
the photon orbits in this last case we have taken into account all the components of the 
velocity of the in-falling plasma, calculated by assuming that it began its in-fall from the 
ISCO. Plunging gas thus preserves its ISCO energy and angular momentum. We found that 
for flares that take place close to the horizon (r^/M ~ 2.1) most photons end up into the 
black hole (as expected); however, the behavior of the response of the escaping photons in 
no case gave us the constant time-lag obtained in the rapidly rotating black hole case, to 
produce QPO features in the PDS. In the case of a Schwarzschild black hole, an (unstable) 
photon circular orbit lies at = 3M < r^s inside the plunging region, from which radius one 
would expect higher-order photons (multiple orbit photons). However, since the plunging 
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X-ray sources also possess significant radial velocity component at this radius close to the 
horizon, many photons emitted at r^/M = 3 are Doppler-beamed in the direction of the 
source motion. The radial beaming effect therefore keeps most of the photons from being 
emitted in an azimuthal direction, effectively suppressing their multiple orbits. 

Finally, to explore even the most likely case to produce QPOs in a Schwarzschild geom- 
etry, we computed the response function of the system for a source at an unstable (Keple- 
rian) circular orbit at r/M = 3.001. At this radius, the high rotational speed of the source 
[v ~ 0.9995c) produces effects similar to the dragging of inertial frames, because most pho- 
tons are beamed in the forward direction by the source rotation. The corresponding response 
function and photon orbits for a source at (pg^i ^ 90° are given in Figures [H^,b, while the 
corresponding PDS and ACF in Figures [9K,b. As it can be seen there, for this specific source 
position, most photons reach the observer after going around the black hole by A0 ~ 37r-|-7r/2 
while the fraction that reaches after A0 ~ n + 7r/2 is comparable but smaller; there are also 
a small number of photons that reach the observer in the counter clockwise direction after 
A(f) ~ 27r + 7r/2; photons that reach the observer after a larger number of rotations are not 
discernible at this resolution. However, as the source phase (ps changes, while the lag between 
the major peaks remains roughly constant, the relative normalization shifts quickly and for 
most of phases the response is dominated by a single peak. This fact, along with the random 
choice of the phase of a given burst of our light curve prescription, "washes-out" potential 
QPO features, leading to the PDS shown in Figure [HK- On the contrary, in the Kerr case, 
the frame-dragging aided azimuthal beaming helps a significant fraction of the photons orbit 
around the black hole for all values of the phase angle 0^. This is an essential difference 
between Schwarzschild and Kerr black hole cases, which we conclude is manifested as the 
QPO feature we see here. 

Up to this point we have only considered the time correlations between the photons 
emitted by the accretion disk without any reference to their energy. However, if these 
photons are of a specific energy (i.e. the Fe transitions observed in accretion disks), there 
exists a relation (redshift) between the emitted and the received photon energies given by 

where Qk is the Keplerian angular frequency, and Vg is the azimuthal component of the X-ray 
source velocity measured in LNRF. It would therefore be of interest to consider, in addition 
to the time correlations, also their additional dependence on the photon energy (or energy 
shift factor) g; the impetus for such a study comes fro m the observed ene rgy dependence of 



the geometry attributed QPOs in GRO J1655-40 (e.g. IStrohmayerll2001al ). We hope to look 



into this as well as the above unresolved issues in a future publication. 

We would like to conclude with some general remarks concerning the "light echo" re- 
sponsible for the QPOs presented in this paper: This is a very generic process and would be 
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present at any source that exhibits a delay in its response function, not necessarily the one 



we described in this note. The same point has been made earlier (iKazanas fc Hua Ill999l ) 
in a different context and we believe that it may have broader applicability. It should be 
noted, in this same context, that the response of an axisym metric disk does not exh ibit any 



such features, however, the response of a warped disk (e.g. Hickox fc Vrtilekl 2005 ) should; 



therefore such disks may exhibit QPO features of the kind described herein at the periods 
comparable to the light crossing time across the disk, features that maybe worth searching 
for in the data. 

Finally, we would like to point out that the lack of apparent phase coherence in this 
model makes the PDS quite noisy and therefore we expect that the detection of such features 
may require their search in the light curves of higher mass objects that would push the 
corresponding frequencies to lower values that are easier to detect (perhaps intermediate mass 
black holes associated with ULXs or nearby AGN), using low background, high throughput 
missions like Constellation-X. 
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5. Appendix 

In this Appendix we show the definition of the local angle 6 of photon emission in 
LNRF and its relation to photon's impact parameter b (which is defined as specific angular 
momentum). We then explicitly derive radial component of the geodesic equation in terms 
of S which we have integrated in producing Figures [2], [Hand [8] in §2 and 3. First, toroidal 
velocity of (either a massive or a massless) particle seen by a local observer in LNRF of Kerr 
metric is given by 



A t 
A 
AV2S 



(10) 

^1 , (11) 



where we restrict ourselves to equatorial trajectories (6* = 7r/2) in the disk plane. In this 
fram e the photon emissi on angle 6 between the propagation direction and the radial direction 
(e.g. iMisner et al.lll973l . p. 675) is defined as 

cot (5 = -r = — -tttaT- ■ (12) 

On the other hand, radial component of the geodesic equation is given by 

E^r^ = R{r) , (13) 

where R{r) = (r^ + — a6)^ — A{b — a)^ for equatorial trajectories and b is photon's impact 
parameter. With the help of equations and (jl]) one can express b in terms of r and S as 

/ f \ 1/2 

b{f,6) = r(^j tan5. (14) 

Substituting equation (fT^ into equation (fT3l) one obtains an explicit expression for r in 
terms of the local angle 6 only as 

pos2 5 (±r A| sec 5| - 2aMAi/2 tan 5) 
""^^^^ " rV2 [r^^r - 2M) + a^{r^ - 2M^) + 2am^ cos(25)] ' ^^^^ 

which is equation ([5]), and the sign in the numerator depends on the direction of photon 
emission in the rest-frame of fiuid/source. Correspondingly, equation ([T^ is now rewritten 

as 

^sin(25)[±rA| sec(5| - 2aMAi/2 tan(5] 
^ ~ 2AV2 [r3(r _ 2M) + a?{r^ - 2M^) + 2a?M^ cos(2(5)] ' ^ ^ 
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Fig. 6. — Timing analysis of random X-ray bursts with a/M — 0.99 for Ts/M — 1.3 (top 
rows), 1.45 (middle rows) and 2 (bottom rows). We set — 6000. Left columns show the 
average power density spectra (PDS) while right columns show the corresponding autocor- 
relation functions (ACF). In the ACF we show two cases: Ni, — 6000 (dark curves) and 100 
(gray curves) for comparison. 
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Fig. 7. — Same as Figure[6]but for a/M = 0.995 (top columns) and 0.9428 (bottom columns) 
with rg = Tms(fl)- 
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Fig. 8. — (a) The response function for a source at an unstable (Keplerian) circular orbit at 
r^/M = 3.001 at a phase 0^ = 90°. The two large peaks correspond to the prograde orbits 
that cover angles = vr + 7r/2 and Stt + 7r/2 respectively. The small peak at midpoint 
corresponds to the retrograde orbit that covers an angle A0 = 27r + 7r/2. (b) The photon 
orbits corresponding to the features of the figure on the left. 




Fig. 9. — Predicted PDS of the light curve of a source at r^/M = 3.001 obtained using the 
response function of Figure [H] and the prescription given in the text, (b) The ACF of the 
same light curve for N\, = 6000 (dark curves) and Nb = 100 (light curves). 



